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In a Peccei-Quinn extension of supergravity the h — > 77 detection rate can be significantly en- 
hanced due to the reduction of the total Higgs decay width. To assess the viability of various 
Peccei-Quinn extensions of minimal supergravity we perform a Bayesian analysis on three such sce- 
narios. The main constraints on these models come from the currently observed Higgs boson like 
state by the Large Hadron Collider and from the WMAP observation of dark matter abundance. 
Our comparative study reveals that under these constraints the PQ violating scenarios with axino 
dark matter are clearly preferred over the minimal supegravity model with the lightest neutralino 
as a dark matter candidate. 
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I. INTRODUCTION 

The ATLAS [1] and CMS [2] collaborations observed 
an excess in the diphoton invariant mass distribution, 
around 126 GeV, with more than 5tr statistical signifi- 
cance. This excess points to a particle with properties 
close to the Standard Model Higgs boson. The observed 
production times decay rates of the particle, for exam- 
ple, are mostly consistent with those of a standard Higgs. 
The diphoton rate, however, is an exception: it is roughly 
twice of the standard value. Although the uncertainties 
are quite sizable, this anomaly suggest a non-standard 
Higgs-like particle. This implies that new physics is af- 
fecting the properties of the newly discovered particle. 

Many new physics models were proposed in the re- 
cent literature to explain the difference of the diphoton 
rate from the standard one. These models assumed ex- 
tra dimensions [3J, a fourth fermion generation j^j, ex- 
tra vector-like leptons 0, a dilation [BJ, extra gauge 
bosons 7J, or multiple Higgs states [8]. Solutions based 
on hybrid models has also been suggested to explain the 
LHC data in Refs. 0. While supersymmetry (SUSY) [TD] 
is one of the robust new physics candidates, ironically, its 
minimal versions are struggling to explain the diphoton 
excess. 

In the Minimal Supersymmetric Standard Model 
(MSSM), for example, it is hard to double the diphoton 
rate of the lightest Higgs boson. The most promising sce- 
nario relies on the existence of light tau sleptons. Within 
the MSSM it is even hard achieve a 126 GeV lightest 
Higgs boson naturally. With top squark masses below 
TeV the lightest MSSM Higgs mass remains below about 
120 GeV. Heavier stop masses, in turn, stretch the hi- 
erarchy between the electroweak and SUSY scales. The 
problem is severe since every single GeV of loop contribu- 



tion to the Higgs mass beyond about 120 GeV requires 
top squark masses further and further above the weak 
scale. These problems have been examined in the context 
of the MSSM [1 1 and its simplest extensions, the next- 
to- minimal MSSM (T5], the natural SUSY scenario [13], 
and the Peccei-Quinn (PQ) extended MSSM [HQS]. 

In a recent work by Blaum et al [TJ] pointed out that 
in the models with a slight PQ violation it is possible 
to raise the Higgs to diphoton decay rates significantly 
depending on the size of the PQ violating couplings in 
the Higgs potential. It is also expected that these mod- 
els can serve better from the point of view of cosmology 
by providing an alternative to the traditional neutralino 
dark matter and leptogensis. With this motivation, in the 
current work, we analyse the PQMSSM in greater detail 
using a Bayesian framework and compare our findings 
with the MSSM to show that the PQMSSM scenarios fit 
the collider and cosmological data with more flexibility. 

The organisation of this paper is as follows. Section [Tl| 
provides a brief overview of the PQMSSM and discusses 
its advantages over the MSSM. Later in the same sec- 
tion also discusses implications on collider phenomenol- 
ogy and cosmology. Section III contains our detailed 
Bayesian analysis of various PQMSSM scenarios. Our 



results and conclusions are given in Sections |IV| and \V\ 
respectively 



II. PECCEI-QUINN VIOLATION IN THE 
MSSM 



To solve the strong CP problem Peccei and Quinn 
(PQ) extended the SM with a global U(l) symmetry 
[IB] [17]. The U(\)pq symmetry is spontaneously broken 
at a scale Apq and the pseudo-Goldstone boson induced 
by this breaking is the axion. The axion mass is related 
to the symmetry breaking scale as 
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The PQ extension of supersymmetric models requires the 
addition of a chiral superfield 



$a = 



V2 
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where s is the scalar axion or saxion, a is the the pseudo- 
scalar axion, a is the axino, and F a is an auxiliary field. 
The scalar and pseudo-scalar axion fields are even while 
the axino field is odd under the R-parity. The masses of 
these field depend on the supersymmetry breaking mech- 
anism. In most cases the saxion is ultra heavy with a 
mass of about Apq while the axino mass is highly model 
dependent. For the supegravity inspired model, which is 
subject of the current paper, the axino mass takes the 
following form|18j. 
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FIG. 1: Ratio of Total decay width and various Branching ratios 
of the SM-like Higgs boson in the PQ-violating mSuGra to the 
conserved mSuGra case for mo = 1000, m 1( i 2 = 500, Ao = (all 
in GeV units), tan/3 = 30, sgn([i) = +. 



Thus depending on the value of k and Apq, m a can take 
values from a few eV to several TeV. This wide range 
of axion masses leads to interesting cosmological conse- 
quences which will be discussed in Section fll B| 

A broken PQ symmetry contributes to the neutron 
electric-dipole moment (nEDM) at tree level [T71 IT9] . 
The current experimental limit on the nEDM is 



d n < 1 1.9 x 1CT 26 | 
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at 90% CL. This limit translates into a lower bound on 
the PQ breaking scale. A model dependent upper bound 
has also been obtained for the PQSuGra case in Ref. |2"0] , 
leading to 



1 x 10 9 GeV < Apq < 5 x 10 11 GeV. 



(6) 
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A. Implications for the Higgs Sector 

With a broken PQ symmetry the MSSM Higgs La- 
grangian can be extended by the following term 



FIG. 2: Various Higgs-Boson observables, R gg xY 
pp^ h 1 susy with respect to ran at the LHC. 
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where | A5 1 < 2. As observed by the authors of Ref. [13] , 
the above term modifies couplings of the Higgs-Boson to 
various SM particles by, 
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From Eqs.( 8apc ) it is clear that the coupling of the SM- 
like Higgs to fermions with T3 = +1/2 receives a posi- 
tive correction due to PQ violation, while for T3 = —1/2 
fermions and the weak-bosons this contribution is nega- 
tive. 
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FIG. 3: Various Higgs-Boson observables, R gg xY 
pp^ h 1 susy with respect to A 5 at the LHC. 



B. Cosmological Implications 

Due to the bounds on the PQ breaking scale in Eq.(|6|, 
the axion mass is always restricted in the range between 
~lx 10~ 5 eV and ~6x 10~ 3 eV. The saxion is ultra- 
heavy and hence is less interesting for our purposes. The 
axino mass can take values between 2 GeV and 1 TeV for 
the case k = 2, and ~ 1 eV - ~ 7 keV for k = 3. Thus 
the axion is always the lightest of the three and hence 
serves as a good hot dark matter candidate. The axino, 
covering a wide range of mass between a few eV to about 
a TeV, can be lighter, degenerate or heavier than the 
lightest neutralino. This makes the PQ violating SuGra 
scenario very intriguing since the both the axino and the 
neutralino can contribute to the cold matter abundance. 
For k = 2, for example, three interesting PQSuGra sce- 
narios are possible depending on the lightest neutralino 
mass, m^ and the axino mass, ma- These are 



PQ-1 (neutralino LSP): m^o < ma, 
PQ-2 (axino-neutarlino co-LSPs): r 



m x'l ' 



PQ-3 (axino LSP): 



m a < m^o. 



In the first (third) case the lightest neutralino xl (axino 
a) is the dark matter candidate. In the second case they 
are both CDM candidates. 

For k > 2 the axino is always the lightest superpartner. 
However, it was found that the only case that is cosmo- 
logically viable for k > 2 is the one with k = 3. Scenarios 
with k > 3 are unable to generate sufficient dark matter 
relic density and hence are less interesting. Therefore in 
the current work we will only analyse the PQSuGra sce- 
narios with k — 2 (denoted by PQ-1, PQ-2, PQ-3) and 
k = 3 (PQ-3'). 



Following Ref . [2UJ [3T] , we calculate the relic densities 
for the axions and axinos by the following formulae 
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Here, f2 Q /i 2 is the axion relic density, Q,^ TP h 2 is the relic 
abundance of non-thermally produced axinos from neu- 
tralino decay and £lj p h 2 is the relic abundance of ther- 
mal produced axinos. 

Axino dark matter lends the PQSuGra model consid- 
erably more viability compared to the minimal SuGra 
model. The properties of axino dark matter, such as its 
abundance and couplings to standard matter, are gov- 
erned by its mass and Apq which are independent from 
the mSuGra parameters. Thus, one expects more flexi- 
bility from a model with more parameters. Obtaining a 
Higgs mass of about 126 GeV, a neutralino abundance 
of 0.22 pc and low fine tuning is impossible in mSuGra. 
The PQSuGra model not only accommodates these re- 
quirements much easier but its Bayesian evidence sug- 
gests that, despite of its extra parameters, it is more 
viable. 



III. NUMERICAL ANALYSIS OF THE 
PQSUGRA MODEL 

In the rest of this paper we consider a mSuGra model 
extended with the PQ violating Higgs coupling Eq.([7|. 
The former is parametrized by the usual four parameters 
and a sign, while the latter adds two more to our full set 
of parameters: 

V = {m ,m 1/2 , A , tan/3, sgn(n),X 5 , A PQ } . (11) 

We use the latest version of SUSY-HIT [H| to calcu- 
late sparticle masses and decay rates, and Mi cr Omegas 
2.4.5 23 to calculate the relic density of lightest neu- 
tralinos. To incorporate the effect of PQ violation we 
modified the relevant parts of SUSY-HIT. 

Due to the extra PQ violating term the Higgs couplings 
to various SM particles receive significant contribution 
compared to the PQ conserving mSuGra case. This, 
in turn, modifies the Higgs branching ratios. To show 
this, in Figure [l] we plot the deviation of the various 
Higgs branching ratios and the total decay width from 
the mSuGra ones. We show this deviation as the function 
of A5 for fixed values of {mo, m-y/2, tan/3, sgn(/j,)} = 
{1000 GeV, 500 GeV, GeV, 30, +1}. 

As expected the deviation rises with increasing A5 
and becomes maximal at A5 = 2. Due to this the 
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total Higgs decay width, for example, becomes about 
10(18)% smaller for A5 = 1(2) compared to the mSuGra 
case (which is equivalent to A5 — 0). This happens 
because for A5 > the h — b — b coupling becomes 
smaller. After folding in changes of other Higgs cou- 
plings the h — > 77 branching ratio becomes larger by 
about 6(18)% for A5 = l(2)compared to mSuGra. There 
are similar changes to the h — > VV* , where V = W ± , Z 
branching ratios. The h — > gg branching ratio increases 
by about 8(19)% for A 5 = 1(2). This translate into 
about (100 - (100 + 8) x (100 - 10)) = 2.8%, and, 
(100 - (100 + 18) x (100 - 18)) = 3.2% decrease in the 
h — > gg partial decay width for A5 = 1, and, 2 respec- 
tively. Since Higgs is dominantly produced through the 
gluon fusion, the ratio 



QBay 
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T mSuGra 
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also decreases accordingly for Higgs production at the 
LHC. It is also clear from the plot that, for the selected 
parameter point, an overall increase of the diphoton event 
rates by about 1.03 — 1.14 is predicted by PQSuGra over 
mSuGra. 

As we saw, in the PQSuGra scenario it is possible to 
enhance the diphoton production rate and dark matter 
is less constraining than in mSuGra. In the light of the 
LHC excess and the strict WMAP dark matter abun- 
dance constraint we aim to quantitatively compare the 
feasibility of PQSuGra and mSuGra. In what follows we 
calculate Bayesian evidences for both models and evalu- 
ate their odds compared to each other. Beyond the LHC 
Higgs search and WMAP we also use data from LEP, the 
Tevatron, and various low-energy experiments. 

The calculation of evidences involves a full scan over 
the parameter spaces of both models. Motivated by nat- 
uralness, this scan is done over the following parameter 
ranges: 

• mo € [10, 2000] GeV 

• mi/2 G [10, 2000] GeV 

• Aq € [-3000,4000] GeV 

• tan/3 e [0,62] 

• A 5 € [0, 2] 

• Ap Q e [l x 10 9 ,5 x 10 11 ] GeV. 

A. Bayes Factor Calculation 



Evidence against the base model 



— 0.5 Not worth more than a bare mention 
0.5 — 1 Substantial 
1-2 Strong 



> 2 



Decisive 



TABLE I: Interpretation of Bayes factors on Jeffery's Scale. 
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The likelihood function can then be constructed as 

1 
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Since theoretical predictions Of" depend on the theoret- 
ical parameters listed in Eq.Q, the likelihood Ci carries 
the same dependence. 

In terms of the likelihood the posterior probability dis- 
tribution for parameter pj is then given by 



f V({x})£ J] dx 



(15) 



where C = TT >Ci- An a priori probability distribution 

i 

■p({x}) is assumed for each parameters. While the PQ- 
SuGra parameters are listed in Eq.(|6]), the mSuGra pa- 
rameter space is smaller: 



V = {■m ,m 1 / 2 , A ,tan/3, sgn(/j,)} 



(16) 



The evidence £ for each model is calculated by inte- 
grating the posterior density over all input parameters. 
To give a specific meaning to the evidence, the Bayes 
factor is constructed: 



QBayes — 1°£ 



10 



£pQSuGra 
£mSuGra 



(17) 



This evidence ratio quantifies the odds of PQSuGra 
against mSuGra. Odds are interpreted in terms of the 
Jeffreys scale as shown in Table [T] 



To assess the relative viability of PQSuGra compared 
to mSuGra we calculate their relative Bayes factor which 
is the ratio of their evidences. First, for each observables 
Oi we calculate \\ using the predicted and experimen- 
tally measured central values, Of 1 ' and £)? xp ', and the 



B. Experimental Constraints 

Our likelihood function includes the following experi- 
mental and observational bounds on various observables 
and the (s)particle masses. 
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1, Precision observables (POs) 
• dp = 0.0008 ± 0.0017 [H], 



• a 



SUSY 



(3.353 ±8.24) x 10" 9 [25 



BR(b -> sj) = (3.55 ± 0.26 ± 5%(th.)) x 10~ 



2. Bounds on sparticle and Higgs boson mass from 
LEP-2/Tevatron 



m^o > M z /2 GeV 2HEl 



In order to estimate various LHC observables we calcu- 
late the SM and SUSY Higgs decay width and branching 
ratios at NNLO + NNLL, wherever available. Numeri- 
cal values for other observables (except for relic-density, 
^dm) were also calculated using SUSY-HIT. For relic den- 
sity predictions we use MicrOmegas. 

To safely avoid a charged LSP we assume that masses 
for all the sleptons and squarks are at least 10 GeV above 
the axino and/or neutralino LSP masses. 

Before closing this section we would like to add that 
for the numerical implementation of various PQSuGra 
scenarios, we assume the mass difference between the co- 
LSPs is within 10%. 



• m x ± > 103.5 GeV [27]. 



IV. 



RESULTS 



• mj > 99 GeV 

• m h > 114.4 GeV [28]. 

3. LHC data on yfs — 8 TeV Higgs searches 
4- LHC Higgs search 
We use the following LHC Higgs search observables 
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Here 7t gg - ry ,7i- gg 2i2v, and, 7t gg u are ratios of diphoton, 
2l2v (where I = e ± ,/i ), and 41 production rates in PQ- 
SuGra relative to the Standard Model. We did not in- 
clude TZ ggT ± T T because there is insufficient data corre- 
sponding to this observable. 

The following are the experimental values used, given 
in the Ref. 0, 

• m h = 126 ± 0.5657(ea;p.) ± 2(th.) GeV 



• R ggil = 1.8 ±0.5 



• Rgg2l2v = 1.4 ±0.5 

• R gg4l = 1.4 ±0.6 

5. Dark matter abundance 



n DM = 0.1123 ± 0.0035 ± 10%(th.) 



The y/s = 8 TcV LHC Higgs decay observables show 
a significant deviation from the SM prediction particu- 
larly for the diphoton, di-lepton and the four lepton cases 
where the observed central value for the former turns out 
to be about 1.5-1.8 times larger than the SM. We have 
shown that PQSuGra can change these observables sig- 
nificantly in the right direction compared to the mSuGra. 
To gain detailed insight in Figs. [2] we plot R ggl7 , R gg bb, 
R gg vv, an d R gg vv with respect to the Higgs mass for 
PQSuGra. These plots are the result of a scan over the 
PQSuGra parameter ranges given above. The plots tell 
that for most of the cases the diphoton event rate drops 
considerably with increasing Higgs mass. Yet there still 
seems to be a narrow region where the Higgs mass and 
Rggn agrees with the data within 2a. (We assume about 
1 GeV theoretical uncertainty in the Higgs-Boson mass 
calculation.) The other two observables, R gg 2\2v and 
R gg u ' are also in good agreement with the LHC measure- 
ments. To show how sensitive each of these observables 
are for the PQ violating effects, in Figs. [3] we plot them 
with respect to A5. This plot tells that it is possible to 
satisfy the LHC data for a wide range of the A5 . 

We calculate posterior probabilities for the SuGra 
and PQSuGra models for two popular choices of pri- 
ors: (a) the flat (uniform) prior which is constant in 
a finite parameter region; and (b) the log prior cx 
(morai/ 2 ) -1 - These posterior probability distributions 
are then marginalized to the various parameters and plot- 
ted in Figs. [4] and [5] As these figures show a large value 
of m > 1.5 and m 1 / 2 > 0.5 TeV, A < -3 TeV and 
tan j3 ~ 20 is preferred by the experimental data. 

The integral of the posterior probability distributions 
allows us to calculate the Bayes factors as presented in 
Table |H| This table along with Table [!] shows that the 
PQSuGra scenario is somewhat better then the mSuGra 
when only Higgs and POs are imposed. However, scenar- 
ios where the axino is dark matter are strongly preferred 
over mSuGra. This is because in the latter it is hard to 
satisfy the LHC Higgs mass constraint, WMAP and g~2 
simoultaneously. In contrast, the PQ — 3 (and PQ — 3') 
models satisfy WMAP in a wider range of the PQSuGra 
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FIG. 4: Posterior probability distributions using the flat prior. 
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FIG. 5: Posterior probability distributions using the natural prior. 

parameter space according to Figs. [3| 



V. 



CONCLUSIONS 



We have preformed a Bayesian analysis of the mini- 
mal and Peccei-Quinn violating supergravity models. We 
compared the viability of the two models in light of the 
LHC Higgs searches at </s = 8 TcV, the WMAP data 
on the relic density of dark matter of the Universe, along 
with data from various other experiments. Our study 
reveals that PQSuGra scenarios with an axino LSP are 



clearly preferred by the collider and astrophysical data. 
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FIG. 6: Posterior probability distributions for the relic density 
(!1dj/) for various PQSuGra scenarios. 



Observables 


Q Bayes 


PQ-l 


PQ-2 


PQ-3 


PQ-3 


Higgs Searches at the LHC 


0.244 


0.305 


0.663 


0.506 


+ POs & LEP data 


0.238 


0.322 


0.724 


0.626 


+ WMAP data 


0.181 


0.231 


1.693 


2.466 



TABLE II: Bayes factors for various PQSuGra scenarios for 

m G [10,2000], m 1/2 G [10,2000], A G [-3000,4000] 
(all in GeV units), tan/3 G [0,62], A 5 G [0,2], \ PQ G 
[1 x 10 9 ,5 x 10 11 ] (in GeV). 
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